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1 Introduction 



Considering quantum mechanical indistinguishable particles, one uses to speak about their 
"statistical" (or "exchange") interaction. It has no classical analogue, since there is no 
classical interaction force, but it influences the many-particle spectrum. Under usual cir- 
cumstances, with bosons and fermions only, the rules of combining single-particle states 
for indistinguishable particles are affected, while particles of different species remain in- 
dependent. For anyons Jl|, [|, however, this interaction is nontrivial, since it leads to the 
absence of single-particle states, and even more, to the possibility that particles of differ- 
ent species can become mutually "entangled"; hence the problem of multispecies anyons 
demands special consideration. 

It is also believed that this problem may be of direct physical interest, particularly in 
connection with the t — J model, which might be relevant to the high-T c superconductivity. 
The question about statistics of excitations therein, spinons and holons, is still unsettled, 
and it is anticipated that their "mutual statistics" may be nontrivial ||. Another possible 
application is the layered Hall media where FQHE quasiparticles in different layers can be 
viewed as anyons with mutual statistical interactions Q . Finally, the same situation arises 
in the usual quantum Hall liquid: In a state where exactly m Landau levels are filled, the 
excitations are characterized by "charge vectors" of dimension m and there are different 
nontrivial phase factors associated with interchanging different sorts of excitations ||. 

In this paper we address certain points concerning thermodynamics of multispecies 
anyons, by starting from microscopic considerations on the wave functions. For the sake 
of simplicity, we do the calculations for two species only, but a generalization to more 
than two species would be straightforward. The fugacity expansion of the thermodynamic 
potential now involves two fugacities, and the equation of state, respectively, two densities; 
hence, the cluster and virial coefficients are labeled by two indices. The mixed virial 
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coefficient of the lowest order, an, is trivially calculable, but exact expressions for higher- 
order ones are certainly not easier to obtain than in the one-species case. To proceed 
further, we investigate the case when all the particles are in the lowest Landau level 
of a strong magnetic field, and, using microscopic arguments as in Q, we obtain the 
equation of state which coincides with the one conjectured in Q and derived in j|] from 
a definition of fractional statistics a la Haldane || (exclusion statistics), using a state- 
counting argument. A mean-field interpretation of this equation is also proposed. Also, 
in this (and only in this) case it becomes possible to interpret the system in terms of 
statistical mechanics of a gas of free particles, by following essentially the same lines as in 
|10(|. A statistical distribution is derived, which again coincides with the one for exclusion 



statistics ||, 11]. We summarize the expressions for the cluster and virial coefficients of a 
multispecies system, in the Appendix. 



2 Mutual fractional exchange statistics 

When particles are distinguishable, their wave functions have no definite behavior with 
respect to their interchange, but one can impose certain conditions that they should satisfy 
under a double interchange, or, equivalently, a winding of one particle around the other 
one. In three dimensions, as usually, wave functions are single- valued, whereas in two 
dimensions a non-trivial phase factor may arise. In fact, if there are M species of particles, 
one has an M x M matrix a a b of statistical parameters, such that the wave function of 
the whole system picks up a phase 

• exp[z7ra aa ] under an interchange of two particles of species a 

• exp [2iira a b] under winding a particle of species a around a particle of species b 
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(provided, in both cases, that the complete closed path encloses no other particles). The 
second condition, of course, holds for a = b as well, but the first one in this case is stronger. 
Note that there is periodicity in a Q a's with period 2, but in a a fe's with period 1. Note also 
that since "winding particle a around particle 6" and "winding particle b around particle 
a" is one and the same operation, in the system of their center of mass, the matrix a ab 
must be symmetric. 

The above condition on the wave function can also be understood within the approach 
where statistics appear as kinematic properties of the configuration space [|J . The config- 
uration space for a system containing M species of identical particles can be constructed 
in analogy with the case of a single species J12[| . Irreducible representations of the funda- 
mental group of the configuration space determine properties of the wave function under 



exchange of particles. In [12| it was found that the fundamental group for the configuration 
space for a system with M species of identical particles is Bj^ 1 ___^ M , a generalization of the 
braid group B^. One-dimensional irreducible representations of BjVi...AT_m (corresponding 
to spinless particles) are labeled by M parameters associated with exchange of particles of 
the same species (which may be identified with the above a aa ) and ^M(M — 1) parameters 
associated with exchange of particles of different species (which may be identified with the 
above a ab ) |||. 

One important observation is now in order. A physical picture usually associated 
with anyons is the model of charge- flux composites ||, where the origin of the anyonic 
interchange phase is the Aharonov-Bohm interaction of the charges and fluxes ascribed to 
the particles. In the multispecies case, each species a would be characterized by a charge 
e a and a flux <f> a . Now, if one writes the Hamiltonian in the regular gauge, including the 
relevant vector potentials, it turns out that a singular gauge transformation to the free 
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Hamiltonian (anyon gauge) is possible only if the equality 

e a 4>b = e-b<t>a (1) 

holds, for each a and b. In this case, one will have 

Oiab = e. a 4> h /2n, (2) 

and the matrix a a b is indeed symmetric. Equations (|l|) and (|2|) together imply that 

a a b = ±^a aa abb- (3) 

Thus, only the diagonal statistical parameters are independent, save for the sign. However, 
a situation may be imagined when there are several Aharonov-Bohm gauge fields. Eq. ([!]) 
then becomes 

44 = 44, (4) 

(3 labeling the gauge fields, while (|2|) becomes a a b = 44 Z^, and (||) no longer holds. 
Therefore, one may consider all a a bS as independent, but in the Aharonov-Bohm model 
this would, generally speaking, demand to have several gauge fields. Note at the same 
time that in a model where the fluxes appear as a result of coupling the particles to 
Chern-Simons gauge fields O], (Q) is satisfied automatically. 

Consider now the simplest case of M = 2, and denote for brevity a\\ = ai, «22 = «2> 
a 12 = 7. A many-particle wave function of a system of Ni particles of species 1 and N2 
particles of species 2 that satisfies the above-discussed interchange conditions reads in the 
anyon gauge 

* = Il(z j -z k ) a * n(c--Cnrn(^-c m ) 7 ^({^} 5 {u}) (5) 

j<k m<n j } m 

01,02,7 > is meant, while the replacement (zj — z^) 01 — ► (Zj — z^)~ ai is to be made 
if ai < 0, etc. $ is a single-valued function symmetric with respect to the coordinates 
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{zj} (species 1) as well as with respect to the coordinates {Cm} (species 2), but arbitrary 
with respect to z <-> £• This is the starting point for constructing the (JVj. + iV2)-particle 
spectrum. Obviously, as far as 7 is not an integer, all (JV1+JV2) coordinates are "entangled" 
together, while if 7 is integer, the problem can be factorized into the one of N± particles 
and the one of N2 particles. 

Knowing the spectrum, one may find the partition function jv 2 = ex P ( — /3-£^Vi7V 2 ) , 
then proceed to the grand partition function 

S= z^Z2 2 Z NlN2 , (6) 
Ni,N 2 

where z a = exp(/3^ a ) are the fugacities, in order to obtain the cluster expansion 

ln~= £ (7) 

fci,fc 2 

and the virial expansion for the equation of state 

ki,k 2 

(see Appendix). 

Clearly, for the virial coefficients a n o and ao n one has 

a n o = a n (ai), a 0n = a n (a 2 ), (9) 

where a n (a) is the n-th virial coefficient of anyons with statistical parameter a. Among 
these, only 0,2 is known exactly 

A 2 , (10) 

with A = yj2irl3/m. The mixed second-order coefficient, an, involves the N\ = N2 = 1 
case. Eq. (B) then turns into 

M/ = (z-C) 7 ^(^C), (11) 
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a 2 (a) 



where <3? is any single-valued function. Since it may always be represented as a sum of 
symmetric and antisymmetric functions, the partition function Z\\ can be written as 

Z n = Z 2 ( 7 ) + Z 2 (l + 7 ), (12) 

Ziici) being the 2-anyon partition function with statistical parameter 7. Using the formu- 
las from the Appendix, we get 

an = a 2 ( 7 )+a 2 (l + 7) =7(1-7)A 2 . (13) 

For 7 = 0, 1, this coefficient vanishes, reflecting the fact that particles of different species 
are mutually independent. 



3 Quantum and statistical mechanics in the lowest Landau 
level 

Calculating higher-order virial coefficients would demand solving a multiparticle problem, 
which is in general impossible. In the one-species case, however, -when the system is 
projected in the lowest Landau level (LLL) of a strong external magnetic field, the whole 
virial expansion can be obtained |J]. Note that for the gap above the ground state to be 
bounded from below and thus the projection onto the LLL to be justified, the particle 
fluxes must be antiparallel to the external field, for example if B < 0, one should have 
4> > 0. If the particle charge e is taken to be positive, then the necessary condition is 
a £ [0, 1] 0. The equation of state then is 

pp = PL ln(l + — —) (14) 
V 1 — av ) 

4 Note that this convention for signs is opposite to Q and coincides with 
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(pL = muj c /Ti is the Landau degeneracy per unit area, p = N/V is the density, v = p/pl 
is the filling factor, uj c = \eB\/2m is half the cyclotron frequency; "strong magnetic field" 
means (3uj c 3> 1). 



of fractional statistics ||, usually referred to as exclusion statistics, with a being the 
quantity by which the number of available single-particle states diminishes as one particle 
is added to the system. Why the equations coincide can be understood within a mean-field 
approximation (where it becomes possible to speak about single-particle states): Since the 
anyon fluxes are antiparallel to the external field, the flux <j) = 2ira/e of a particle added 
partially screens the external flux <3? = BV and thus diminishes the number of states 
available for other particles; when N particles are added, the Landau level degeneracy per 
unit area becomes pl — pa. 

A multispecies generalization of fll4| ) within the framework of exclusion statistics is 
as follows: If a ao is the decrease in the number of single-particle states available for a-th 
species caused by addition of one particle of 6-th species, then 



with v a = pal PL % %■ 

The physical picture would be the following: a particle of a-th species carries statistical 
charges e£, which couple to the statistical gauge fields, and an electric charge e which 
couples to the external magnetic field and has to be common to all species. In fact, we 
are forced to restrict ourselves to the Landau density pi being the same for all species 
(whereas in ||], where mean field arguments are used, different Landau densities pLa 
are possible), because it is only in this case that exact eigenstates can be constructed. 

In the two-species case, by using microscopic arguments for the LLL spectrum, we will 



It has been shown H that the same equation (14) is obtained from Haldane's definition 




(15) 
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now show that the above equation indeed holds with a a bS being the mutual statistical 
parameters defined previously: it is then almost certain that (|i~5| ) is correct for any number 
of species. 

If there is only one statistical charge, then, as already emphasized, 7 is not independent 
from the aw's. Now, both fa and fa have to be positive, and it follows then from (|l|) that 
e\ and e 2 must have the same sign; choosing it again to be positive, we have a%, a 2 G [0, 1], 
and, by virtue of (|2|), 7 is positive as well, so the constraint @ would read 

7 = \/aia 2 . (16) 

However, our considerations will hold in the case of several gauge fields as well, and the 
result will be true for arbitrary 7 G [0, 1] (with the condition, however, that the gap above 
the ground state be bounded from below by a quantity of the order the Landau gap 2to c , 
which is certainly true for 7 small enough). 

By analogy with ||, we add a harmonic attraction X/^i"^ 2 ^muj 2 ^ and observe, from 
(||), that the spectrum of the (iVi + A^-particle system in the LLL becomes that of a 
system of N\ and N 2 bosons (mutually independent) in a "harmonic-Landau" potential 
whose single-particle spectrum reads 

si = uj t + (u t - lo c )£ (17) 



(u>t = y/ujf+uj 1 , I = 0, 1, . . .), plus a constant shift 



A£ , Ar 1 AT 2 (ai,a 2 ,7) 



NtiNt-l) N 2 (N 2 -1) ' 
■a\ -I a 2 + NiN 2 ^ 



(tut - u c ) (18) 



2 2 
coming from the anyonic prefactor in (||). One then has 

Z Nl N 2 = e-P AE ^^> a ^Z b Nl Z b N2 , (19) 
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where 



Zn (1 _ e -/3(ut-Wc)\h - e -2/3(ut-u>o)\ . . . (1 - e -NP(wt-f> c )\ ^ 

is the iV-boson partition function in the harmonic-Landau potential. 

For the cluster coefficients, one gets (see Appendix for the prescription for passing to 
the thermodynamic limit) 

*« = '^n(i-^) ■«-*■*■ ■ p« 

/i=i '2=1 

(22) 

and bok 2 = bk 2 o\a 1 ^a 2 - Note that, as expected, (^) does not reduce to (|2l|) when &2 = 
(in particular, 6^1 fc 2 certainly vanishes for 7 = 0, while b^o is 7 independent). The virial 
coefficients calculated herefrom take the form (/3u c 3> 1) 

Ofafa = ~ ^./^^"^' {[K " l)" 1 " «M7 fe2 + [(«2 - l) fc2 - «(H7 fc1 } . ( 2 3) 



and the equation of state then is 

PP = PL 



ln(l + - +ln 1 + 



(24) 



l — — 7^2/ V 1 — CX2V2 — i v \. 

which indeed coincides with (Jl5|) . 

A mean-field interpretation is possible in this case as well. Indeed, what figures in the 
Hamiltonian for a-th species is the sum eA + J2p e a^> where A is the vector potential of 
the magnetic field and A@ is the /3-th gauge field potential created by all other particles. 
If all the gauge fields are averaged, so that B@ = J2b Nb4>b/V = J2b Pb4>\,-> then the Landau 
level degeneracy per unit area becomes 5j(|ei3| — X)/3 e a^) = PL — J2bPb a ab, which is 
precisely what leads to (|i~5|). 
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As well as for the one-species case, this equation is valid only when v\ and vi ar e small 
enough — in fact, both denominators have to be positive; one of them vanishing means 
that the critical filling is achieved for the corresponding species. In the special case fllq) , 
the condition is 



a\v\ + V«2^2 < , ; =?■ (25) 
\/max(ai, a<i) 

4 Interpretation in terms of single-state distributions 

We now wish to interpret the exactly solvable model of anyons in the LLL in terms of a gas 
of free particles with a particular statistical distribution. We will use arguments similar 
to those of [pi] for the Calogero model. Description of single-state partition functions for 



multispecies systems follows Ref. [11 



Let there be single-particle states labeled with i, common for all species a, but with 

(i) 

energies e a possibly depending on a. The particles in different single-particle states are 
assumed to be statistically independent, which implies that the grand partition function 
is a product of single-state grand partition functions 

s=n s(4) 

i 

with HW depending only on the Gibbs factors x a = e^ a ~ E ^: = SW({x a }). The 
statistical distributions (average number of particles of different species in state i) are 

n»=z a ^-hiS» (27) 

We assume that In may be expanded in integer powers of the Gibbs factors (again, we 
specify the formulas containing the series expansions for two species): 

oo 

lnS»= £ f kl kA'4 2 - (28) 

fci,fc 2 =0 
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It follows then from (I 



4°= E KfkrkA 1 ^ 2 - (29) 
fci,fc 2 =i 



Combining ( |2"6| ) with (^) , summing over single-particle states and comparing with (^) , 
we get the relation 

kfcifc 2 = fk 1 k 2 cr k 1 k 2 i (30) 

where 

^ lfc2 = E e_/3(fcl£ " )+fc24l)) - (si) 

i 

We use the many-particle partition functions for anyons in the LLL in a harmonic well 
(|lS|), with ui — ► 0, to calculate the coefficients fk 1 k 2 i n the expansion p8|). Expressing the 
cluster coefficients in terms of the many-particle partition functions (but not yet using 
the thermodynamic limit prescription given in the Appendix), taking into account the 
expression ( |l7| ) for the single-particle spectrum, in order to calculate ak 1 k 2 i an d then 
using (|3l|) and (|30|), we eventually obtain 

fk l0 = ± nV-^), (32) 
f klk2 = -7^-^ n . a — T — ) n a — — )■ (33) 

(l = l 12 = 1 

We thus conclude that in the limit u — > 0, multispecies anyons in the LLL can be 
described by the partition functions (p8j) and the statistical distributions ( p9| ) with the 
coefficients of their expansions in powers of the Gibbs factors given by (|32])-(|33|). 

Now we consider directly the case u> = 0. All the single-particle states have then the 
same energy u c , therefore all the relevant quantities do not depend on the index i any 
longer, and summation over single-particle states reduces to multiplication by the number 
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of states in the LLL, piV, e. g. InH = pi,V lnH^ 1 ). It turns out then that H^ 1 ) can be 
represented as 

^=Y[E a (34) 

a 

with 



2„ = 1 + 



(35) 



1 - J2b a abnb 

so that for the coefficients of the expansions of the single-state grand partition functions 



fcl,fc 2 



we obtain 



*W = n 1 



Zi=2 



fei 



Pi ;fclfe2 = -7^ 11(1 ) 11 ( 1 " - 



h=2 



k 2 -l 

n 



(36) 

(37) 
(38) 



and P2\k\ki — -fljfefci |o!i<-+a2 ■ 

Eqs. (|32| ) and (|37]) have been previously derived for fractional statistics in one di- 
mension, starting from the Calogero model (see also a recent derivation directly for 
exclusion statistics [15 1). Eqs. (33) and ( |3"8| ) are generalizations onto the multispecies case. 

The partition functions H a , in addition, satisfy the equations 



2 tt = l + a:«n s 6 



S„.b— ot n 



Substituting (35) into (p9|), we get the equation for the statistical distributions n a 



n a /x a = 1J(1 - a bc n c ) aab {l - n a - ^ 



ab c n c 



:$ab — a ab 



(39) 



(40) 



These statistical distributions coincide with the ones for exclusion statistics ||, |ll]] 
after the identification 

9ab = a ab (41) 
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(in our model, consequently, the exclusion statistics parameters are symmetric, g a b = Qba)- 



It is now easy to calculate the cluster coefficients, using (30). Using the thermodynamic 
limit prescription at each order k\ + &2, ( |3l|) becomes 



Cfcife — PL 



y e -f}(k 1 +k2)uic 



(42) 



and, with the use of (^)-(33), the formulas ([H|)-(|2"2|) are recovered. Note finally that 
once (|34j) and (^) have been correctly "guessed" , the equation of state is nothing but the 
thermodynamic identity (3PV = InS. 



5 Concluding remarks 

We have studied quantum and statistical mechanics for mutual fractional exchange statis- 
tics. We have presented a microscopic solvable model, anyons in the LLL, which, on 
one hand, involves mutual fractional exchange statistics and, on the other hand, repro- 
duces the statistical mechanics for mutual fractional exclusion statistics. This establishes 
correspondence between the two approaches to mutual statistical interactions. 

We hope that due to the similarity between anyons in the LLL and FQHE quasipart- 
icles, the present work will help a better understanding of excitations in quantum Hall 
systems. 

In addition, intimate links between anyons in the LLL and the Calogero model suggest 
that the microscopic solutions for multispecies anyons found in this paper can be used to 
construct solutions for a multispecies Calogero model which would reveal mutual exclusion 
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Appendix 

We summarize here, for reference purposes, the formulas concerning the cluster and 
virial expansions for a two-species gas, and give the relations between Zn 1 n 2 , b ni n 2 , and 
a nin2 , up to the fifth order. 

The grand canonical partition function is 

3 = E ^^ 2 Z NlN2 (A.l) 

N lt N 2 

and the cluster coefficients bk ± k 2 are determined in terms of Z^^'s by writing the loga- 



rithm of ( A.l ) in the form 

ln~ = E (A-2) 

ki,k 2 

The pressure and the densities are 

respectively (V is the area); writing 

/W= £ "mA'p* 2 , (A.4) 

ki,k 2 

using ( |A.3j ) with In 5 substituted from ( |A.2| ) , and matching coefficients at equal powers of 
z's on both sides, one gets the equations which allow to express a^k 2 i n terms of bk^s. 
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The explicit formulas are given below. For the cluster coefficients, one has 
ho = Zw] 



ho = - 2^10 + ^20, 
hi = —ZqxZ\q + Z\x\ 



ho — 3^10 ~~ ^10-^20 + ^30; 

hi = Zq\Z\q — Z\qZu — Z01Z20 + Z21; 



ho — —^Z 10 + Z 10 Z2o — -Z 20 — Zi Z 30 + Z 4:0 , 

hi = — Zqi Zf + Z\qZ\\ + 2^01^10^20 — Z11Z20 — Z10Z21 — Z01Z30 + Z31, 

*^ 2 2 2 ^2 2 

^22 = —-^Z 01 Z l0 + Zq2Z 10 + 2ZqiZiqZh — -Z u — ZiqZi 2 + Z 0l Z2Q 

— Z02Z20 — Zq\Z2\ + Z22; 



ho 
hi 

hi 



1 



^10 ~~ ^10^20 + Z\qZ 2 q + Z 10 Zso — Z20Z30 — Z10Z40 + Z50, 

v3 v o v v2 v 1 o v v v 1 v v2 1 v2 



5 

ZqiZiq — Z'i Zu — 3ZqiZ{ Z2o + 2^10-^11-^20 + -^01^20 + ^10^21 



— Z20Z21 + 2Zq\Z\qZ^q — Z11Z30 — ZiqZ^i — Z01Z40 + Z41, 

o /^2 v3 v *73 o v rj2 v 1 v v2 i v2 v o v2 v v 
^ Z 01 Z 10 ~~ Z 02^10 ~~ -J^Ol^lO^H + z 10 z ll + ^10 Z 12 — -jZ 01 ZioZ20 

+ 2Z02Z10Z20 + 2Z01Z11Z20 — Z12Z20 + 2Z01Z10Z21 — Z\\Z 2 \ 

— Z\qZ22 + Zoi-^30 — Z02Z30 — Z01Z31 + Z32. 

If a harmonic oscillator regularization is used, the thermodynamic limit is meant as u> — > 0, 



with a particular prescription at each order of the cluster expansion [17]. A treatment 
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analogous to the one given in [18] shows that to get the correct thermodynamic limit for 



bk 1 k 2 , the necessary prescription is to identify (f Cl +k^)/3mui' 2 
the virial coefficients, defining 

bk!k 2 = &fcifa/&io&oi) 
a>k x k 2 — a k 1 k 2 / v j 

one has 



V, in the limit uj — > 0. For 

(A.5) 
(A.6) 



aio 



CL20 

an 



-ho, 
-hi; 



030 = 4&20 - 26 30 , 

an = bji + 4bnho - 2b 2 i; 



a 40 = -20&20 + 18^20^30 - 36 40 , 

% = - 66^620 - 24611^0 + 3611621 + 12620^21 + 9611630-3631, 

022 = -96o26n - 3bfi + 6611612 - 12602611620 - 9611620 + 6612620 + 6602621 
+ 6611621 - 3622; 



a 50 = 1126^ + 32640620 + 186| - 1446 30 6l - 46 50 , 

041 = 6fi + 86?i6 20 + 4Gb 2 n b 2 20 + 1606n&l " ^ 2 u hi - 24b u b 20 b 2 i - 806l 6 21 
+ 26|i - 126^630 - I2O611620630 + 246 2 i6 30 + 46n63i + 246 20 6 3 i 
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+ 16&H&40 - 4&41, 

a 32 = 16602611 + 66^ -126^612 + 48602^1^20 + 32^620-32611612620 

+ 646o26ii6^ + 806^6^ - 32&i 2 &lo - 246 02 &ii&2i - 206^621 + 86i 2 6 2 i 

- 32602620621 - 64611620621 + 86^ + 8611622 + I6620622 - 24602611630 

- 246^630 + 12612630 + 8602631 + 12611631 - 46 32 ; 

with the obvious symmetry with respect to permutation of the subscripts. 
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